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Abstract 

In this paper, we consider self-adjoint difference equations of the form 

-A(a n _iAy„_i)+fc„y„ = Xy n ,n = 0, 1, . . . (0.1) 

where a n -\ > for all n > and b n are real and X is complex. Under the as- 
sumption that a n - 1 satisfies certain growth conditions and is limit point (that 
is, the associated Hamburger moment problem is determined), we prove that 
the existence of an exponentially bounded solution of (10.11) implies a bound 
on the distance from X to the spectrum of the associated self-adjoint opera- 
tor, and that if a solution of (10.11 ) is bounded by a power of n for n sufficiently 
large, then X G o(B). Here, B is a certain self-adjoint operator generated by 
(10- 1 b - These results are the difference equation version of differential op- 
erator results of Shnol'. We use this to then prove that the spectrum of 
the associated orthogonal polynomials contains the closure of the set of X 
for which we can find a polynomially bounded solution. We also present a 
result concerning the invariance of the essential spectrum under weak per- 
turbation of the coefficients. 

Keywords: difference equation; orthogonal polynomials; exponential bound; 
spectral theoiy 



1 



INTRODUCTION 



1. Introduction 

In this paper we consider the second order difference equation in self-adjoint 
form 

-A(a n -\Ay n -y) + bny n = Xy n ,n = 0,1,... , (1.1) 

where a„_ i > 0, and b n is real for n = 0, 1 , . . . , and A is a possibly complex 
number. Here and throughout this paper, A is the forward difference operator 
defined by 

4?« =gn+l ~gn- 

We shall also assume that for any solution y n of (11.11 ). y_i =0. Equation 
(11.11 ) may be written in the form 

- a„y n+ \ -a n -\y n -y + (b n + a n + a n -\)y n = Xy n 

and by making the change of dependent variable y n h-» (— l) n y n , we may 
write 

a n y n +i +a„-iy n -i + (b n + a n + a n -\)y n = Xy n (1.2) 

Using a theorem commonly attributed to Favard (see Ismail [6], p. 31), 
we see that there is a sequence of orthonormal polynomials {p(n; X} which 
satisfy (11.2b and there is a measure dfl(X) such that 




p(n;X)p(m;X)dn(X) = 8, 



Thus a solution of (ll.lb is then (— l) n p(n; X). Note that the polynomials 
{p(n; X)} satisfy the initial value problem consisting of (11.2b and the initial 
conditions p{— 1; X) = and p(0; X) = 1. There have been a number of pa- 
pers in the past few years devoted to the study of orthogonal polynomials, 
particularly the interaction between the recurrence coefficients in (11.2b . the 
orthogonality measure d/J.(X), and the polynomials {p(n; X)}. For recent 
references, see J6]. 

There is an intimate connection between the spectral theory for (ll.lb . or- 
thogonal polynomials, continued fractions, asymptotic expansions of Stielt- 
jes integrals, moment problems, variational inequalities, oscillation and dis- 
conjugacy, and a host of other topics The three volumes by Henrici [4 ] are 
perhaps the best self-contained reference for the connections between or- 
thogonal polynomials, continued fractions, Pade approximation, and some 
of the other topics mentioned above. Disconjugacy and oscillation for (ll.lb 
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is covered in Kelly and Peterson ifTTTl . In an earlier paper ifloll . we discussed 
the connection between the oscillation of solutions of (11.11) and the essen- 
tial spectrum of a linear operator generated by (ll.lt ; for the definition of 
this linear operator, see below. For other results on the oscillation of so- 
lutions of (11.11 ) and orthogonal polynomials, see van Doom [19] and [20]. 
Mingarelli lfT2l discusses second order difference equations from the point 
of view Volterra Stieltjes integral equations, and unifies many of the results 
on oscillation, spectrum, and comparison of solutions proved separately for 
(11.11 ) and second-order ordinary differential equations in formal self-adjoint 
form. Other general references involving (ll.lt are Ismail [0 on orthog- 
onal polynomials and related topics, and Jones and Thron [|9] on continued 
fractions and their relationship with Pade approximation, moment problems, 
and orthogonal polynomials. We have certainly not included all references 
to books or papers which contain something on these topics; that would be 
a rather large task. 

In this paper, we pursue the connection between the spectral theory of (ll.lt 
and exponential and polynomial upper bounds on solutions of (1.1). It is 
important to note that we do not require estimates of the orthogonal polyno- 
mials, but that some solution of ( 1 . 1 ) have an appropriate estimate. We prove 
the difference equation versions of some theorems of Shnol' and Simon; see 
Glazman [3], p. 175-182 for proofs of Shnol's results in English in the ellip- 
tic differential equation case and lfl5l for Simon's results. A special case of 
these results were announced in section 2 of [16 | and [18], but the proofs of 
the results stated here are given here for the first time. The results given in 
lfl6l and lfl"8l are very restrictive, since the author had to assume that (1.1) is 
nonoscillatory, a„_ i = 1 for all n, and b n is bounded below. In this paper, we 
remove all the restrictions we just mentioned, but place growth constraints 
on the a n . 

We shall need the relevant operator theory, which is found in Hinton and 
Lewis @. Let H = f{Z+) be the Hilbert space of square-summable se- 
quences with inner product 

oo 

n=Q 

The associated norm will be denoted by || . ||. Let the linear operator B :H i-> 
H be defined by 

D(B) = {y£H:By£H} 
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where 

(By) n = -A(a n -\Ay n -\) + b„y„ 

The resolvent set of B is the set 

p(B) = {X G C : B - X G M(H)}, 

(B(H) is the set of bounded linear operators on H) and the spectrum of B is 
the set g(B) = C \ p (B). An eigenvalue of B is a (necessarily) real number 
X such that the equation By = Xy has a non-zero solution in H. An eigen- 
value is called a discrete eigenvalue if the associated eigenspace has finite 
dimension (that is, the eigenvalue has finite geometric multiplicity), and the 
set of discrete eigenvalues of B is denoted by G c i(B). The essential spectrum 
of B is the set o ess (B) = o(B) \ Od(B). Using the spectral theorem, we may 
characterize these parts of the spectrum in an alternative way: the essential 
spectrum is the set of real numbers which are of infinite geometric multi- 
plicity or are limit points of the spectrum; the discrete spectrum is the set 
of points of the spectrum which are isolated. For more information on these 
terms and the spectral theory of linear operators, see Weidmann |[2Ti . 

Finally, equation (11.11) is said to be limit-circle if every solution of (11.11) with 
Im(X) 7^ is in H; (11.11 ) is said to be limit-point if it is not limit circle. Note 
that this classification is independent of X (see [ Q). For a discussion of these 
terms and the reason behind the names limit-circle/limit point, see Chapter 
1 of Akhiezer JT]. An important result for us is the following proposition, 
which is a combination of Theorem 2 in Hinton and Lewis [5 ] and Corollary 
2.2.4 and Theorem 2.1.2 in Akhiezer 0]. 

Proposition 1.1. The linear operator B is self adjoint if and only if (11.11 ) 
is limit point if and only if the associated Hamburger moment problem is 
determined. 

In the limit-circle case (or if the associated moment problem is not deter- 
mined), the situation is more complicated. In that case the operator B is not 
self adjoint by Proposition 1.1, but a restriction of B may have self adjoint 
extensions. Welstead ([22] and [23 ]) and Shi and Sun lfl4l have shown that 
a boundary condition at infinity is needed to construct the self adjoint exten- 
sions of this restriction of B. For the abstract theory of self adjoint extensions 
of symmetric operators, see Weidmann [21]. 

2. The Main Results and Discussion 

The main results are as follows. 
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Theorem 2.1. Suppose there are constants C\,y,riQ,L > such that for all 
n > no 

n 

\Aa„-i\ < Cia„_i, £ a l-\< LeJn 

a«cf suppose d 1 - 1 b z's /imif poi«/. Suppose there are numbers j3 , C2 > aw<i a 
solution y n of (11.11 ) smc/z that for all n > «o 

a«<f A real. There there is a constant C = C(X,C\,C2) such that 

1 /2 

d(X,a(B)) <c(e 2p -lj 

where d(X,o(B)) is the distance between X and the spectrum ofB. 

Theorem 2.2. Suppose that the first two hypotheses ofTheorem \2.1\ are true, 
but suppose there is a solution y n of (11.11 ) such that for any j8 > there is a 
finite constant C3 (/$ ) > with 

W\<CniP)<P n 

for all n> no and that X is real. Then X G o(B). In particular, if there is a 
real number 6 and a constant C4 {X ) such that 

\y n \ <C 4 n 9 

for all n > no, then X G o(B). 

Theorem 2.3. Suppose there are constants C\,y,no,L>0 such that for all 
n > no 

n 

|Aa„_i| <Cia„_i, a 2 k _ x <Le yn 

k=n Q +\ 

and suppose (11.11) is limit-point. Let 

E = {X G R : has a polynomially bounded solution}. 
ThenE c a{B). 
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We must again point out that these theorems do not require an estimate on 
the orthogonal polynomials obtained from (ll.lt . Also, these results should 
be compared with those of Smith lfl6l and lfl"8l . Our results here do not re- 
quire a hypothesis on the b n (other than that it is a real sequence) and the 
hypotheses on the a n are weak enough to allow for many interesting exam- 
ples, such as polynomial ore exponential behavior. 

Theorems 12.11 12.21 and 12.31 should be compared to the differential equation 
results found in (starting on p. 176) and [151 - In [3], the potential q of 
an elliptic differential operator —A + q is required to be bounded below, but 
this condition is weakened in Ifl5l . A key part of the proof of Theorem l2.2l is 
Theorem 2.1 (see below), which gives an a-priori bound on A;y„_i in terms 
of y n . There are corresponding theorems for the differential equation case 
in 0, P- 176-178 and |[T5l . The proof of Theorem I2.3l relies on a general- 
ization of Theorem 9, p. 197 of and is related to Stieltjes conjecture (see 
section 58 of 0). Several years ago, Simon [15] proved a theorem like 12.31 
for a class operators related to the N-body problem of quantum mechanics. 

The reader may suspect that the assumption of self-adjointness in Theorems 
12. 1 [ 12.21 and I2.3l is an artifact of our method, and may be removed by using 
other methods. But self-adjointness seems to be essential to prove theorems 
such as these, even in the differential equations case. The reason for this is 
as follows. In the limit-circle case it is known that self-adjoint extensions of 
a restriction Bo of the operator B may be considered as finite-dimensional 
perturbations of Bo and thus (by Weyl's theorem) have the same essential 
spectrum (see [3], p. 10, 44-46). But discrete eigenvalues of an abstract 
operator may in fact disappear if the operator is perturbed by an arbitrarily 
small bounded operator (see Kato iflOl ; the theorem is known as the Weyl- 
Von Neumann Theorem). Indeed, one of the earliest quantum mechanical 
models exhibits the phenomenon of loss of discrete eigenvalues - the Stark 
effect (a hydrogen atom in a uniform electric field; here the perturbation 
is unbounded, see Simon). Thus, if discrete eigenvalues are present, the 
spectra of the various self adjoint extensions of a restriction of B may be 
quite different. Because of this it is not clear whether the results of this 
paper hold in the limit-circle case. A condition we do not investigate is the 
essential self-adjointness of B. These results would probably require us to 
know some form of relative boundedness for (ll.lt . We do not know if such 
results yet exist in the literature. For a discussion of relative boundedness 
and its connection with essential self-adjointness, see Weidmann |[2T1 
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3. Proof of the Main Results 

The following theorem plays a key role in the proof of Theorem 12.11 as does 
the corresponding result for the differential equations case; see Simon |fT51 . 

Theorem 3.1. Suppose that there is a constant C\ > such that for all n 
sufficiently large 



|Aa n _i| < C\a n -i 



(3.1) 



and suppose y n is any sequence ( not necessarily a solution of (11.11 ) ). Then 
there is an integer m, and for integers r, s, and n which satisfy m<r < s <n, 
we have 



£at 1 (A W _ 1 ) 2 <2(l + (C 1 + l) 2 ) £ aUyl 



k=r 



k=m—l 



Proof. First, choose m large enough such that 

| Aa^t— 1 1 < Cifl/t-i for k>m—\. 
For k>m— 1, (13.31 ) implies that 

cik 



(3.2) 



(3.3) 



1 



or 



1-Ci < 



au_ 

ak-i 



< 1+Ci. 



(3.4) 



Now 

fati(Aw-i) 2 <I«ti(lw-il + N) 2 - 

k=m k=r 

By first squaring on the right of this last inequality, and then using the in- 
equality lab <a 2 + b 2 on the cross-terms, we have 

n n— 1 n 

E 4-1 (Ay*-i) 2 < 2 £ 8*^ + 2 

&=m k=m— 1 

In the first some on the right of (|3.51 l use 



(3.5) 



k=m 



cik = a-k-v 



a-k 
au-\ 



and (1341 ) to get 

£aLi(AM-i) 2 <2(l + (C 1 + l) 2 ) £ fl 2 y 2. (3.6) 



k=m 



k=m—l 



□ 
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We would like to note here that the proof of 13. H is much shorter and easier 
than the proof given in [18] and the proof for the differential equation case 
on p. 176-178 in Also, we do not require y„ to be a solution, as does O, 

ed, and na. 



3.1. Proof of Theorem I27T 

The proof of Theorem 12. 1 1 relies on using [J/,V]z = U(Vz) — V(Uz), the 
commutator of U and V, which is defined for z £ D(U)nD(V). In the theory 
of the Schrodinger equation, commutator estimates such as those used below 
have been found to be very useful; see Simon [15]. The results are in Section 
C.4 of 1151 . Commutator estimates have also been used in the theory of 
Jacobi matrices; see (H. The proof of Shnol's theorem found in Simon iflBI 
is where we learned the idea of using the commutator. The complication 
we have with the difference equation case in this paper is mainly due to the 
product rule, which is 

A (s„t n ) = s n+ iM n + t n As n 



or 

A (s n t n ) = s n At n + t n+i As n 

We also define the shift operator (Eu) n = w„_i and the norm on the integers 
in the interval [s,t] by 

t 

II l\2 V 2 

n=s 

Proof of Theorem \2J] Let hq be a positive integer chosen such that for all 
n > no + 1 

|Aa„_i| < Cia n _i and \y n \ < C 2 e P ". 
Let r — 1 > no and let v r be the sequence 



V r — V r n 



1 if no < n < r, 
otherwise. 



In what follows, A is differencing with respect to n. Now v r y is in the domain 
of B (because it has finite support), so 

[B,v r ] (v r y) = B (v r y) - v r By 



or 



B(v r y) = [B,v r ]y + Xv r y. 
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by using the fact that y satisfies (B — k)y = 0. Thus, 

(B-X){v r y) = [B,v r }y 
Using the first product rule above, we have 

[B,v r ]y n = -A (a„- 1 A (v n „_ \y n ))+ b n v r , n y n + (b n -X)v m y n . 
Expanding this out by using the product rule again gives us 

[B , v r ] y„ = - (v r>n A (a„_ i Ay n - i ) + a n Ay n Av r . n + Av r . M A (Cl n — iyn— 1 , 

+ tfn-lJn-lA V r „_i 



+ (^-^)v ri „y„. (3.7) 

But in this last equation, the first and last terms on the right add together to 
give zero, since y is a solution of (ll.ll ). Thus 

[B,v r ]y„ = -Av r> „ (a n Ay n + A(a B _iy„_i)) - a n ^y n ^ x A 2 v r ^ n 

Using the product rule again for A gives 

[B,v r ] y n = -Av m (a n Ay n + y n Aa n _ \ + a„_ i Ay„_ \)-a\_ t A 2 v r ,«- 1 • 

We need ([B, v r ]y n ) 2 . This is an easy, but tedious operation. We first square 
both sides, and in the cross terms, use the inequality lab <a 2 + b 2 to get 

([B,v r ]y n ) 2 < C (Av r , n ) 2 \a\ [(Ay n ) 2 + {Ay n ^f\ +y 2 n _ l (Aa n ^) 2 } 

+ Cfltiii(AVi) 2 - (3-8) 

Now we may write this using the norm as 

(\B,Vr}y„) 2 < c{\\aAy\\ 2 rr+3] + \\a(AEy) ||f r , +3] 

+ \\A(Ea)Ey\\l r+3] + \\EaEy\\ 2 r _ lr+3] } (3.9) 

Note that we have used the following facts 

supp (A (Ev r )) C integers in the interval [r — 1 , r + 2] ; 
supp (A 2 {Evrf) C integers in the interval [r— 1, r + 3]; 
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in obtaining (13.81) from (13.91 ) and the constant C > in these last inequalities 
does not depend on n. These may be derived from the definition of the 
sequence v r . Now, by changing the index of summation, 

h A y\\lr+3] = \\EaAEy\\l +l r+4] < \\EaAEy\\ 2 [rr+4] . (3.10) 

Also, by the hypothesis on a„_i, 

HaAEyllf^ < (d + l) 2 \\EaAEy\\ 2 rr+3] (3.11) 

Using first (13.101 ) and (13.1 II ) in (I3.9I ). and then using Theorem j:\vice, we have 

|| [B,v r ] y\\ 2 < K { \\Eay\\f r _ lr+4] +2\\Eay\\l r+3] + \\Eay\\j r _ u+3] } , 

(3.12) 

where Kis a positive constant. But (13.81) allows us to write this last equation 
immediately above as 

II (B - X) v r y\\ 2 < K { ||£ay||f r _ 1/+4] + 2||^ll[, r+3 ] + ll^llfr-i.r+3]} • 

(3.13) 

We now define, for r > no + 1, 

F(r)= t al_ x y\. (3.14) 

n=«0+l 

Using (I3.7I ). (13.121 ) may be written as 

\\(B-X)(v r y)\\ 2 <K(4F{r + 4)-3F(r-2)) (3.15) 

or 

|| (B - X) (v,-y) f < K { (F (r + 4) —F(r — 2)) (3.16) 

Using the exponential bound on y„ gives us that for any 8 > there is an 
ro > «o + 1 such that for r > ro 

F(r)<eW +S > £ a 2 n _ v (3.17) 

We now show that for a given 8\ > there is a sequence of integers r p — > oo 
such that 

F(r p + 4) <e 2 P +Si F(r p -2) (3.18) 

(strict inequality is important here as we shall see below). If no such se- 
quence exists, then there is an r\ > no + 1 such that for all r > r\ 

F(r + 4)>e 2p+S[ F(r-2) (3.19) 
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From this it follows that there is a constant M > such that for all r large 
enough 

F(r) >Me {2p+Sl)r . (3.20) 
Using (13.181 ) and (I3.22I ). and the hypothesis on Y, a n-v we nave 

e {8 x -8-i)r < _L for r > maxr0iri . (3.21) 

We now choose 8i > 8 + y, which gives us a contradiction. Thus, our as- 
sertions involving (13.181 ) is true. Now translating this back to the norm and 
operator gives us that 

\\(B-X)w p \\ 2 <K l {e 2 ^-\), (3.22) 

where w p = w pn = v r jtl y n . Now w p converges weakly to zero but ||w p || = 1, 
so w p is not compact but is bounded. Using Theorem 10, p. 14 of Glazman 
01 on (13.221 ) proves the theorem since 8 > is arbitrary. □ 

Proof of Theorem |Z2] The first statement follows directly from Theorem 
12.11 To prove the second statement, we need only show there is an expo- 
nential bound, and apply the first statement of the theorem. Now for any 
P > 0, if n > no, then 

\y n \<C^=C^ 

But this last term is bounded independently in n. □ 
The following theorem is proved in Glazman J3], p. 197, Theorem 9. 

Theorem 3.2. For any £ > and for jj,-almost all values of X (where fl is 
the spectral measure ofB), there is a constant L(e, X) > such that for all n 

\(-l) n p(n;X)\ <L{e,X)n l / 2+e . 

Here the p(n;X) are the orthogonal polynomials obtained from (I1.2I ). 

But we need more than this to prove our conjecture that under the hypothe- 
ses of Theorem 12. II E = o(B). 

To show these results are non-trivial, we consider a special case of an exam- 
ple found in Wimp 1241 . Consider the difference equation 

-n(n+\)y n+ i+2n 2 y n -n(n-l)y n -i = -nXy n . (3.23) 
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In ll24ll it is shown there are two linearly independent solutions of (13.231 ) 
such that 

y hn = n-^e 1 ^ i {\+o{\)), (3.24) 
J2,n=^ 3/ V 2 ^ 1 (l+ (1)) (3.25) 

We shall change the dependent variable using the idea of Hinton and Lewis 
l5l . For a difference equation of the form 

— Pnyn+l + qnyn ~ Pn-\yn-\ = kc n y n 

let 

yn = \/~Cnyni 

_ -1/2 -1/2 

b n =c~ l q n . 

Then equation in y is of the form dl .2b . Note that this change of variables is 
related to the concept of equivalent continued fractions; see [9]. Applying 
this to ( |3.23l > gives the equation 

- (n(n+\)V 2 y n+l +2ny n - (nin-l)) 1 / 2 ^ =-X% (3.26) 

which has two linearly independent solutions 

y l ,n=n-^ 4 e 2 ^ X (\+o(l)), 
y 2 ,n=n- l l A e- 2 ^{l+o{\)) 

The sequence a n = (n(n + l)) 1 / 2 satisfies the hypotheses of Theorem 12.11 
and (13.261 ) is limit-point (see 1 1 ]). For all A there is a bounded solution. 
Our results show that the spectrum of the associated self-adjoint operator is 
(—00,00). This result appears to be new. 



4. Invariance of the Essential Spectrum 

In this section, we consider the difference equation 

- A [(«„_! + TJ„_i) Ay n -l] + On + Wn)y n = (4. 1) 

as a perturbation of (ll.U . Here, we require that a„_i +f] n -\ > for all 
n. The result we present is applicable even if (11.11 or (I4.lt do not generate 
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self-adjount operators or whether the associated moment problems are de- 
termined or not. 

To state the theorem, we need to extend our operator-theoretic discussion 
found in the introduction using ||5]. Let A : H G H deonte the operator de- 
fined using ( 14.11 ), that is, 

D(A) = {y£H:Ay£H} (4.2) 

where Ay n is the left-hand side of (14- 1 b - We now define the minimal opera- 
tors Aq and So as follows. Let T = A or B (the same discussion works for 
both) and let 

D\ = {y £ H : only finitely many of the y„ are non-zero}. 

Now let 




be the restriction of T to D\. A calculation shows that T\ is symmetric and 
densely defined, so it has a closure Tq. This minimal operator may or may 
not have self-adjoint extensions. For the relevant terms and theorems which 
we have used here, see Weidmann ETI . 

Theorem 4.1. Suppose b n > a > Ofor all n sufficiently large and 

lim^ = lim^=0 

«^°° a n n^oo Jj n 

This result is the difference equation analogue of the result for differential 
equations found in Glazman Q- The proof of Theorem 14. 1 1 is not given 
here, since it closely resembles the proof of Theorem 25 in section 8 of l3l. 
The key to the proof is to show that Ao — Bo is relatively compact with re- 
spect to Bo and then to use Weyl's theorem on the invariance of the essential 
spectrum under relatively compact perturbations; see Weidmann [21] and 
Chapter 1 of 0. 
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